1-30-1998 11:34AM FROM GATECH ECE/CSIP <.348948363 

' 

E2I X73 : FlNAl REPORT 


P.2 


N A £) - / — / t f 


Genetic Adaptive Control for PZT Actuators 

Jeongwook Kim, Shelley K. Stover and Vijay K. Madisetii y sS ~ 3> ~ O 

/S/ts/?/*/? — 206816 ^ 7'*// 


Abstract — A piezoelectric transducer (PZT) U capable of 
providing linear motion if controlled correctly and could 
provide a replacement for traditional heavy and large servo 
systems using motors. This paper focuses on a genetic 
model reference adaptive control technique (GM3AC) for a 
PZT which is moving a mirror where the goal is to keep the 
mirror velocity constant. Genetic Algorithms (GAs) arc an 
Integral part of the GMRAC technique acting as the search 
engine for an optimal PID controller. 

Two methods are suggested to control the actuator in this 
research. The first one Is to change the PID parameters 
and the other is to add an additional reference input in the 
system. The simulation results of these two methods are 
compared. 

Simulated Annealing (SA) is also used to solve the prob- 
lem. Simulation results of GAs and $A are compared af- 
ter simulation. GAs show the best result according to the 
simulation results. The entire model is designed using the 
Math works' Simulink tool. 



Keywords — GAs, GMRAC, PID, PZT, Simulink 


Figure 1. PZT System 


I. Introduction 

Most linear actuators are comprised of motors which are 
very heavy, but are capable of providing stable linear mo- 
tion over long periods of time. Although the PZT actuator 
can provide a highly controllable linear motion, the inher- 
ent characteristics of the material can change due to aging 
and temperature. Therefore, an adaptive control technique 
is desirable to compensate for these changes. The GMRAC 
technique which incorporates genetic algorithms was used 
to set the PID parameters and thus control the velocity 
movement of a mirror being driven by the PZT. 

Genetic Algorithms (GAs) have been useful in solving 
various kinds of problems. GAs are primarily used to solve 
optimization problems which have a complex, nonlinear 
search space. The design automation for controllers is one 
application for GAs. There arc several Instances where 
GAs have boon applied to this type of problem: Lee and 
Takagi designed a fuzzy system using GAs for the inverted 
pendulum [5]; Michalewicz used GAs for solving the opti- 
mal control problem [7]; and Das, Goldberg, Maclay and 
Dorey used GAs for system identification [1,6]. 

In this research, two different optimization approaches 
are proposed. The first is to change the PID gain param- 
eters and the other is to adjust an additional input value. 
The performance of the two different methods are com- 
pared and simulated annealing (S A) is also used to compare 

Thi* nn cwch wu supported in pen by NASA L«ngl«y - Flight Elec- 
tronic* Division (1VV4-1995). Submi**ion on 1 Oct IMS. 

Tlit Author# *re with the DSP I.*bor*tory, Electric tl ond Computer 
Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0350, 
USA. vkm9ec-gMcch.edu. httpi//www.ee fetech.edu/u#ere/2lS. 


the performance with GAs. According to the simulation 
results, the GA outperformed the other methods. Specif- : 
ically, the first approach using GAs for PID gain control ; 
had the best result. 

II. PZT System 

The PZT system described here consists of stacks of 
piezoelectric structures which were invented at NASA's 
Langley Research Center. When an input voltage is applied 
to the PZT, a force is exerted by the piezoelectric which 
delivers a displacement to the attached mirror. These de- 
vices can be planed mechanically in series and connected 
electrically in parallel. The motion of each element adds in 
phase to produce the total motion which is almost equiv- 
alent to the sum of motion of each element. The system 
configuration is shown in Figure 1 [10]. The real photo 
image for this system is in Figure 2. 

In this simulation, the PZT plant system was modeled 
in Simulink. The entire system model is shown in Figure 3. 

III. GAs and GMRAC 

According to evolution theory, nature will adapt to a 
changing environment in an efficient manner. The infor- 
mation about each individual is contained in its genes, and 
the processes of natural selection and progressive breed- 
ing build a strong population [2, 3]. Natural selection en- 
sures that the genes from a strong individual are present 
in greater numbers in the next generation than those from 
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Figure 2. Photo Image for PZT System 


&SL 




Verity 




no 




PZTdrhw 


1 


1&2M 






rdx*y amatd 
U ) cottmdl 


r* j elctiiyawj 

s 


veheilyorar 


" H ebdiy | 




late 


L H fora 




Acafanra 


Figure 3. PZT Controller 


a weak individual. After a number of generations, this 
process can combine stronger genes to produce an even 
stronger individual whose genes will strengthen the popu- 
lation. The Genetic Algorithms (GAs) are based on this 
natural selection and genetics. They combine survival of 
the fittest among 6tring structures with an information ex- 
change to form search algorithms with some of the innova- 
tive flair of human search. In every generation, a new set 
of strings is created from pieces of the old. GAs efficiently 
exploit historical information to determine a new search 
point with improved performance [4]. 

Complex structures can be coded by using simple rep- 
resentations such as bit strings (4]. These act like natural 
genes. A number of these strings or individuals form the 
population of GAs. In the next step, successive genera- 
tions of individuals will be built by using genetic operators. 
These genetic operators are used to make new individuals 


Figure 4. Basic Procedure of GAs 


by combining the genes of other individuals. When this 
process is controlled properly, it can produce a rapid im- 
provement in the overall fitness of the population. This 
produces a population of structures that is well-adapted to 
the problem. After number of generations, the structures 
represent good solutions: The fittest individual should be 
an optimal solution. GAs can perform the optimization 
process even in large and complicated search spaces. The 
basic procedure of GAs is shown in Figure 4. The most 
popular method of representing the problem is a binary bit 
string. In this string, we can use either ‘0’ or '1' as a bit 
value. GAs using this method are known as Binary-Codtd 
Genetic Algorithms (BCGA) [8]. 

Choosing the fitness evaluation for GAs is a very impor- 
tant factor. The fitness function must be devised for each 
problem to be solved. Given a particular string, the fitness 
function returns a single numerical “fitness". For many 
problems, particularly function optimization, it is obvious 
what the fitness function should measure. But this is not 
always the case. It is necessary to calculate the fitness of all 
individuals in an entire population. Therefore, it is better 
to use as simple a function as possible. 

The GMRAC uses GAs with the plant model [9]. Re- 
ceiving plant output and reference input information, GAs 
select the best PID gains at each time step from the candi- 
date populations. It is very important for the adaptation 
to occur between one time step and the next. After a set 
number of time steps, a new output is generated using the 
best set of PID gains from the simulate time steps. There- 
fore, the process may not converge before the output must 
be set, but a near optimum value is selected for best per- 
formance. 

The block diagram for the GMRAC is shown in Figure 5- 
GMRAC needs two types of plants, the real system and 
the plant model are needed for the implementation. In 
the simulation stage, the PZT simulink model represents 
both systems. A disturbance is used to simulate the real 
environment and adds noise to the system. 
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Figure 6- GMRAC Operation Module for P2T Control 


ation and mutation are used as genetic operators, linear 
scaling is used for preventing early convergence, and elite 
passing is used to keep the best string during the entire 
process. 


IV. GA IMPLEMENTATION 

The GMRAC operation module for the PZT system is 
shown in Figure 6. In the simulation, two types of opti- 
mization methods are used. The first one is to find the 
optimal set of PID gains. The other is to find an addi- 
tional reference input to the error signal. The first method 
involves finding four gain factors but the second method 
just needs one value. Therefore, the size of the chromo- 
some string is smaller in the second case. The simulation 
is performed for both methods. The model of the first case 


The detailed block diagram of the GMRAC module for 
the PID gain optimization is shown In Figure 8. In this 
module, the fitness function gives the best fitness values 
to the set of PID gains which have the smallest difference 
between the reference input and the system output. 

The second case Is shown in Figure 9. In this case, the 
additional reference input is added to the reference input 
part in the original system. The size of string is normally 
smaller than when using PID gain control and thus reaches 
the optimal value much faster. The GMRAC module for 
this case Is shown in Figure 10. 


S 
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Figure 11. Simulation Environment 



Figure 8. System for Additional Input 



Figure 10. CMRAC Module for Additional Input 


V. Simulation Result 

The environment for simulation is shown in Figure 11. 
In Figure 11, there is a subblock for the GMRAC routine. 
Two output graphs are shown bottom right side. The sim- 
ulation time and method can be decided by the simulation 
menu from the menu bar in the simulator. 

Four types of fitness functions were used to get the best 
result. The definition of each of the fitness functions are as 
follows: 

a e(Jfc) - Vm (k + 1) - y(k+l), e d (k) = e(*) - c{k - 1) 

• Pi » c<t(k), Pi « e(*) 

a type 1: pa 

a type 2: 50pi + 500pa 

• type 3: 500p2 

a type 4: 50Qpi + 50p2 

• Fitness : eack \„ c 

The objective is to reduce the error. Therefore, the op- 
timization value must have the smallest error. The mini- 
mization problem is easily converted to the maximization 
problem. This minimization problem is changed to the 
maximization problem by using the fitness value as . 
The definition of each case is summarized in Table 1. Cases 
1 to 8 use the PID gsun control. Case numbers which are 
greater than 20 use the additional reference input control. 
Three input sets are used to compare the performance be- 
tween these two methods. The input value starts at “0” 
and reaches “1” at 1 msec in the first case. The first 18 
simulation results use this step input command. Simula- 
tion results with different fitness function types are shown 
in Figure 12. As shown in table 1, the noise is not added 
in these cases. The next simulations are for additional ref- 
erence input control. The same types of fitness functions 
are used in this method, and the results are shown in Fig- 
ure 13. The simulation results are almost the same for both 
methods. The fitness type does not affect the simulation 
result in these cases. 
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fitness 1 

No-noise 


fitness 2 

No-noise 


fitness 3 

No-noise 

Case 4 

fitness 4 

No-noise 


fitness 1 



fitness 2 


Case 7 

fitness 3 

With-noise 

Case 8 

fitness 4 

With-noise 

Case 9 

No GAs 

No-noise 

Case 10 

No GAs 

With-nolse 

Case 21 

fitness 1 

No-noise 

Case 22 

fitness 2 

No-noise 

Case 23 

fitness 3 

No-noise 

Case 24 

fitness 4 

No-noise 

Case 25 

fitness 1 

With-noise 

Case 26 

fitness 2 

With-noise 

Case 27 

fitness 3 

With-noise 

Case 28 

fitness 4 

With-noise 


Tabl* 1 . Definition of Each Cate 



Figure 12. Simulation Result (cate 1-4) 


The output results of GAs and SA with type X fitness 
are shown in Figure 14. GAs are not applied to case 9 
and noise is not added. The simulated Annealing (SA) 
technique was also simulated for comparison using the same 
cases. Simulation results where white noise is added to the 
system are shown in Figure 15. 

The second set of simulation uses a different input set. 
The input velocity command starts at value “1" and keeps 
this value to the end of simulation. In this case, the white 
noise is not used. The noise value starts with “0" at 0 sec 
and is changed to “0.1” at 0.15 sec, “-0.1” at 0.2 sec, and 
finally “0” at 0.35 sec. The algorithm uses the fitness type 
1. The simulation result without noise is shown In Fig- 



Fix ur« IS. Simulation Raault (ca*« 21*24) 



Figure 14- Simulation Result (GAs and SA) 


ure 16. The results with noise arc shown in Figure 17. In 
this case, the type of fitness function affects the simulation 
results slightly. The fitness fonction type 1 provides the 
best result when compared against the others. 

The final simulation uses a step function as an input 
Signal. The result is shown in Figure 18. In the Figure, 
section (A) shows the result for the entire region. Sections 
(B),(C), and (D) show detailed results for specified regions 
in section (A). In this simulation, GAs and SA are used to 
solve the problem and four different results emerge from 
the simulation. The legend in the graph is explained as 
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Figure 15. Simulation Result with NoU« 
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Figure 1ft. Simulation Raault without Nob* 


follows. 

• GA1 - GAs for PID gain control. 

• GA2 - GAs for additional input control. 

• SA1 - SA foT PID gain control. 

• SA2 - SA for additional input control. 

According to the simulation results, GAs and SA show 
better performance than the original PID controller. Specif- 
ically, GAs for optimising PID gain control performs the 
best. 



Tim* (**e) 


Fleur* 17. Simulation Result with Noise 


VI. Conclusions 

GMRAC techniques h&ve been shown to outperform other 
adaptive control methods for the control of a PZT actuator 
system. Good control of PZT systems would enable them 
to be used in place of heavy, noisy motors in a variety of 
applications. 

According to the simulation results ; GAs and SA show 
better performance than a normal PID controller. Fur- 
thermore, GAs for PID gain control show the best result 
yielding a fast and stable response. GAs and a plant model 
were combined to form the GMRAC model, and it demon - 
strates good performance in dynamically controlling the 
system with or without noise. 
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Figure 16- Simulation Result with Four Methods 
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